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PART  I : 


MICROPROCESSES 


I.  PROJECT  SUMMARY 

During  this  last  year  one  project  is  being  completed  and 
two  ma j or  projects  have  been  undertaken. 

Past  observations  with  free-fall  microstructure  recorders 
which  sample  temperature  and  salinity  on  a millimeter  scale 
along  a vertical  path  through  parts  of  the  ocean  are  being 
analyzed.  Extensive  mi.crostiucture  profiles  were  collected 
by  Michael  Gregg  at  several  locations  in  the  Pacific  during 
the  past  2-1/2  years.  Although  he  left  Scripps  in  August 
1974,  he  is  :ontinuing  the  processing  of  this  data  here  with 
a programmer,  Norio  Shioura  . He  has  made  several  visits  to 
supervise  the  work  and  expects  to  complete  the  work  on  approxi- 
mately 70  tapes  by  the  summer  of  1975.  The  analysis  being 
conducted  at  Scripps  will  include  the  plotting,  spectra,  and 
deconvolution  programs  that  in  the  past  have  been  applied  to 
a more  limited  data  set. 

When  completed,  this  work  promises  information  on  the 
central  South  Pacific  and  the  Equatorial  Undercurrent  in  August 
on  the  North  Pacific  thermocline  and  on  the  mixed  layer  in 
several  seasons  of  the  year. 

New  analysis  techniques  for  this  data  inventory  will  be 
developed  by  him  at  the  University  of  Washington. 

The  two  new  projects  are  (1)  development  of  the  MSR-10, 
a free- falling  temperature,  velocity,  electrical  conductivity 
microstructure  recorder,  and  (2)  development  of  ship-lowered 
conductivity,  temperature,  depth  (CTD)  fish  for  documentation 
of  temperature,  salinity  and  density  structures  down  to  a 
resolution  of  10  centimeters.  The  instruments  are  being 
developed  to  study  the  velocity  and  density  structures, 
particularly  at  the  boundaries  of  intrusive  features  in  the 

ocean . 


II.  TECHNICAL  REPORT 
1.  MSR-10 

A six  foot  long,  ten-inch  diameter  aluminum  pressure 
casing  was  built  with  four  wings  of  six  foot  radius,  to  provide 
rotational  and  falling  speed  stability.  This  instrument  is 
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basically  similar  to  previous  ones  but  larger,  to  permit 
greater  data  collection  and  storage  facilities,  and  designed 
for  high  stability  and  constancy  of  falling  speed  (see  Figure 
1),  to  permit  measurement  of  velocity  mic r o s t r uc t ur e in 
addition  to  teuperature  and  salinity. 

In  December  1974,  the  instrument  was  taken  out  to  the 
San  Diego  Trough  for  extensive  engineering  tests.  The 
electronic  package  measured  the  rotation  rate  by  means  of  a 
flux-gate  saturated  core  magnetometer  (with  accuracy  to  better 
than  1/2  a degree  orientation  to  the  earth's  magnetic  field), 
pressure,  and  the  vertical  and  horizontal  components  of  tilt 
and  acceleration  through  two  accelerometers.  Figure  2 shows 
a record  typical  of  the  engineering  drops.  As  can  be  seen 
from  the  record,  the  overall  tilt  of  the  instrument  from  the 
true  vertical  was  typically  not  more  than  0.1  degree.  This 
represents  a departure  of  the  instrument  probes  at  the  bottom 
of  the  instrument  by  less  than  3 mm  from  a plumb  line  to  the 
top  of  the  instrument.  The  results  are  in  agreement  with  a 
numerical  model  which  simulates  the  motion  of  the  MSR  through 
a shear  layer  in  the  ocean. 

In  addition,  vertical  accelerations,  corresponding  to 
changes  in  the  falling  speed,  were  found  to  be  smaller  than 
the  resolution  of  Che  vertical  accelerometer,  or  less  than  a 
change  in  1 mm/sec  over  a 10  second  time  per ' od  . 

These  results  encourage  us  to  believe  that  the  measure- 
ment of  velocity  nicrostructure  will  be  possible  with  a 
sensitivity  capable  of  detecting  turbulence  and  shear  in  the 
very  quiet  ocean  environment.  We  plan  to  go  to  sea  in  March 
of  this  year  with  a modified  electronic  package,  already  built, 
providing  ten  data  channels.  Thi  first  look  will  measure, 
in  addition  to  instrumental  rotation,  fall  rate,  tilts  and 
accelerations,  the  gross  or  overall  temperature  variations, 
and  two  sets  of  velocity  and  temperature  microstructure  probes. 
One  set  will  be  mounted  on  the  nose  to  measure  vertical 
temperature  gradients  and  velocity  structure,  an',  two  mounted 
on  one  of  the  four  wings  to  measure  the  horizontal  temperature 
and  velocity  mic ro s true t ur e . In  this  measurement  the  features 
of  the  1SR-10,  namely  a stable,  vibrationless,  free-fall 
platform,  will  be  essential. 

Each  record  from  a drop  will  contain  a continuous  FM 
recording  of  the  above  parameters  for  a 100  meter  vertical 
section,  with  an  estimated  velocity  sensitivity  of  better  than 
0.1  mm/sec,  and  spatial  resolution  of  about  1 mm.  Since  we 
measure  the  thermal  mi cr o s t r uc t ur e in  the  same  manner  as 
previously,  this  first  look  will  allow  s to  answer  some 
important  questions  regarding  the  velocity  structure  of 
thermally  active  and  quiet  regions. 
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In  addition  to  these  preliminary  recording  techniques 
we  are  constructing  a digital  recording  system  for  the  MSR , 
incorporating  a microprocessor  to  permit  programming  the 
instrument  to  examine  interesting  regions  it  encounters  in 
far  greater  detail  than  is  presently  possible,  in  addition 
to  allowing  many  more  data  channels.  We  expect  to  begin 
sea  trials  with  this  electronic  package  in  the  latter  half 
of  1975. 

I 


I 


i 

I 


2.  Conductivity-Temperature-Pressure  (CTD) 

Instr .'mentation  for  Ship  Work 

We  have  taken  a CTD  of  Neil  Brown's  design  to  sea 
several  times  in  the  past  year,  and  have  obtained  over  90 
'yo-yo'  profiles  of  temperature  inversion  regions  off  the 
San  Diego  coast.  In  each  set  of  profiles  the  instrument 
is  rapidly  raised  and  lowered  through  a vertical  section 
of  about  50  m to  study  the  details  of  the  density  structure. 
Emphasis  has  been  placed  on  studies  of  intrusive  features. 

In  one  study  a fine  structure  temperature  recorder  was 
strapped  to  the  CTD  and  two  channels  of  thermal  microstructure 
obtained  for  the  many  profiles.  Analysis  of  this  data  is 
u'derway,  and  it  appears  that  we  can  infer  density  gradient 
structures  to  a higher  degree  of  sensitivity  than  was  possible 
before.  It  has  been  found  that  the  conductivity  sensor 
(see  Figure  3)  on  the  CTD  has  too  rapid  a time  response  for 
the  0.1  second  sampling  interval,  thus  introducing  unnecessary 
aliasing  and  apparent  spiking  in  density;  a new  probe  meeting 
the  correct  sample  requirements  is  under  construction.  We 
plan  to  add  two  thermal  mic r o s t r u c t ur e channels  to  the  CTD, 
and  have  the  entire  system  interfaced  directly  with  our  NOVA 
computer  system.  The  interface  has  already  been  built. 

This  will  allow  us  to  make  density  computations  and  plots  of 
data  at  sea,  hence,  facilitating  adjustment  of  data  collection 
depths  and  experimental  procedures. 

A joint  investigation  with  J.  Cairns  and  personnel  at  NUC , 
utilizing  a towed  thermistor  chain  towed  at  6 knots  by  one 
ship,  and  our  'yc-yo'  capability  with  the  modified  CTD  on 
another  ship,  has  been  discussed,  and  is  tentatively  planned 
for  this  summer.  This  will  provide  the  first  detailed 
examination  of  an  intrusional  feature,  both  spatially  with  the 
towed  thermistor  array  and  temporally  with  our  CTD,  making 
repeated  profiles  of  a 'tagged'  region,  over  a period  of 
several  days  . 

Two  papers  have  been  completed  and  submitted  to  DEEP  SEA 
RESEARCH  and  the  JOURNAL  OF  FLUID  MECHANICS.  The  first  is  a 
paper  concerning  a laboratory  and  numerical  modelling  study  of 


1 


I 
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HI.  ABSTRACTS 

1.  "Design  Considerations  of  Wing-Stabilized  Free-Fall 
Vehicles",  by  A.C.  Mortensen  and  R.E.  Lange 


This  paper  reports  the  results  of  a design  study  of  a 
free-fall,  s e 1 f -r o t a t ing  oceanographic  instrument  utilizing 
wines  to  control  vertical  speed  anu  enhance  stability. 

Numerical  modelling  shows  this  vehicle  to  be  stable,  exhi  iting 
tilts  of  less  than  1 0 “ 2 radians  under  fall  into  -shear  layer 
moving  at  20  cm/sec.  For  non-stalled  wings  the  fall  rate,  u 
within  possible  choices  of  parameters  consistent  with  a Reynolds 
number  of  ' 50,000  is  found  to  be  expressed  by 


where  to  — rotation  rate 
L = wing  length 
C = wing  chord 
mg  = vehicle  overweight. 


With  a vehicle  overweight  of  5 kilograms 
meters,  and  wing  chord  of  .2  meters  the 
fall  velocity  attainable  is  8.7  cm/sec. 


, wing  length  of  2 
calculated  minimum 


"On  the  Decay  of  Grid  Generated  Turbulence  in 
Stratified  Salt  Water",  by  R.E.  Lar.ge 


wl.  lihb.  suppreosion  by  .be 

and  the  transition  to  a system  of  layers  nndnlatinguith 
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scaled  according  to  an  overall  Froude  number  U/LN,  scaling 
the  vertical  wake  width,  where  U is  the  grid  speed,  L the 
mesh  size  of  the  grid,  and  N the  Va  isa  la-Brunt  frequency  of 
the  fluid.  It  is  found  that  the  wake  reaches  a maximum 
vertical  width  at  about  1/2  a Vaisala's  period  and  scaling 
in  vertical  width  proportional  to  the  length  scale  (UL/N) 
1/2  and  then  generally  diminishes  in  vertical  width  until 
about  three  Vaisala's  periods,  when  the  wake  is  scaled  by 
the  length  scale  (U/N)-  The  behavior  of  the  variance 
spectra  of  salt  fluctuations  is  examined  in  time,  for 
different  overall  Froude  numbers  and  are  found  to  redden 
according  to  the  non-dimensional  decay  time  Nt . Finally 
the  equipar t i t ion  of  energy  between  potential  and  kinetic 
energies  is  examined  in  one  setting  of  U,  L,  N space,  and 
eq ui par t i t ion  is  observed  to  hold  by  1/2  a Vaisala's  period 
corresponding  to  the  maximum  wake  growth  point.  A simple 
wake  model  is  presented,  based  on  the  turbulence  becoming 
suppressed  when  the  wake  reaches  a characteristic  Ozmidov 
length,  and  the  time  and  vertical  width  values  predicted  by 
the  model  compared  to  the  observed  values. 
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PART  i I : 


MIDWATER  THERMAL  STRUCTURE 


I.  PROJECT  SUMMARY 

This  report  supplements  previous  reports  on  the  development 
and  use  of  a freely  floating,  yo-yoing  midwater  instrument 
capsule.  A summary  of  accomplishments  for  the  period  Jily  1, 
1974  - December  31,  1974  is  giver.. 


II.  TECHNICAL  REPORT 

1.  Internal  Waves 


During  June  - July  1974,  internal  wave  records  of  21 
days  duration  were  obtained  at  a depth  of  350  m in  an  area  500 
miles  southwest  of  San  Diego.  The  analysis  of  these  data  is 
nearly  completed.  The  resulting  vertical  displacement  spectra 
decrease  generally  as  w 2 between  the  local  inertial  frequency 
(u)  . ) and  the  Brunt  - Vaisala  frequency  (N)  , and  drop  sharply 
outside  these  limits.  Just  below  N there  is  a pronounced 
spectral  peak.  Vertical  coherences  measured  over  separations 
up  to  100  m are  nearly  independent  of  frequency  between  and 
N.  Within  the  range  0 < < . 5N  coherence  decreases  linearly 

with  vertical  separation.  Phase  spectra  indicate  equal  upward 

and  downward  transport  of  wave  energy.  An  internal  wave  model 
proposed  by  Garrett  and  Munk  is  in  general  agreement  with  these 
observations.  Both  the  displacement  and  coherence  spectra 
suggest  that  internal  wave  energy  is  concentrated  primarily  in 
the  first  few  vertical  modes.  These  results  will  be  submitted 
for  journal  publication  in  the  second  quarter  of  1975. 

No  internal  wave  measurements  are  scheduled  for  the 
immediate  future.  However,  as  3 prelude  to  water  motion  studies 
planned  tentativel)  for  1977,  the  study  of  capsule  dynamics  is 
continuing.  In  early  spring  of  1975  the  behavior  of  the  capsule 
in  an  unstratified  fluid  (Lake  Tahoe)  will  be  examined.  Then 
again  in  fall  of  1975,  its  dynamics  in  a thermally  stratified 
lake  will  be  studied. 

By  understanding  the  capsule  dynamics  we  hope  in  later 
experiments  to  be  able  to  infer  water  motions  from  the  corre- 
sponding motion  of  the  capsule. 
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2 . Ocean  Microstructure 

Analysis  of  the  June  1973  mi c r o s t r uc t ur e data  has  been 
completed.  Results  are  reported  in  Williams  (1975),  the 
abstract  of  which  follows: 


A freely  drifting,  midwater  float  has  been 
developed  and  used  to  obtain  profiles  of 
temperature  gradient  (-  .02  m resolution)  and 

temperature  ('.5  m resolution)  at  '6  min 
intervals.  Twenty  houis  of  data  obtained  at 
550  m depth  500  km  offshore  from  Southern 
California  are  analyzed  for  thermal  structure 
and  internal  wave  motions.  Intense  micro- 
structure activity  is  episodic.  Three  micro- 
structure lenses,  one  of  which  is  bimodal,  are 
found:  two  are  associated  with  small  intrusions 

of  '5  m vertical  extent  and  *100  m horizontal 
extent  in  the  horizontal  direction  of  sensor 
motion.  Tne  third  lens  is  interpreted  as  being 
the  result  ot  a shear  instability.  The  Cox 
number  is  '1.5  orders  of  magnitude  greater  for 
profiles  containing  inversions  than  elsewhere. 

Two  microstructure  lenses  are  associated  wit 
internal  wave  shear  maxima,  the  third  is 
associated  with  a shear  minimum.  This  result 
hints  at  a possible  relation  between  internal 
wave  shear  and  microstructure  activity. 


A journal  article  to  more  widely  disseminate  these 
results  is  near  completion. 

In  early  1975  additional  m ic r o s t r uc t ur e measurements  are 
planned  “ a’deep,  fresh  water  lake,  a setting  in  whrch  salinity 
effects  on  the  density  stratification  are  minimal.  °eean 
measurements  are  planned  for  later  in  the  year  with  salinity 
sensors  augmenting  the  present  sensor  complement.  Th, e added 

sensors  will  reduce  ambiguities  present  in  temperature  data 

alone . 
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Part  III: 


SURFACE  WAVES  AND  NEAR  SURFACE  EFFECTS 


I.  TECHNICAL  REPORT 


Work  has  continued  in  the  analysis  of  air  sea  interaction 
measurements  taken  during  the  ARPA  sponsored  OWAX  cruise  of 
FLIP  in  1973. 

Comparison  of  direct  and  indirect  measurements  of  fluxes 
of  momentum  sensible  heat  and  latent  heat  were  carried  out  by 
Dr.  Gregory  Dreyer  in  his  Ph.D.  dissertation  and  will  be  the 
subject  of  forthcoming  publications.  Major  platform  motion 
corrections  were  carried  out  successfully  for  the  first  time 
so  that  the  effect  of  this  error  on  the  direct  covariance 
estimate  could  be  determined.  It  was  found  that  tilting 
motions  of  FLIP  cause  errors  as  high  as  30%  in  the  direct 
measured  momentum  flux,  but  are  negligible  for  latent  and 
sensible  heat  fluxes. 

Analysis  of  high  frequency  response  temperature  sensor 
signals  taken  during  OWAX  by  Steve  McConnell  show  significant 
departures  from  universal  similarity  spectra  forms  for  hig 
Reynolds  number  scalar  fields  with  corresponding  Prandtl 
number.  A summary  report  was  given  at  the  APS  Fluid 
Dynamics  Division  meeting  in  Pasadena  in  November  1974  . It 

was  found  that  for  unstable  conditions  in  the  marine  boundary 
layer  the  temperature  spectra  exhibit  power  law  benavior  less 
"steep"  than  the  expected  f ” 5 ' 3 inertial  subrange  expected 
sometimes  approaching  ~ f'4/3.  Inertial  subrange  constants 
evaluated  when  f_5/3spectra  appeared  were  quite  variable  and 
generally  large,  as  though  a high  frequency  mode  of  dissipation 
were  present.  McConnell  observed  an  apparent  correlation  with 
absolute  humidity  for  about  eight  different  runs  in  four 
expeditions,  and  suggested  a possible  explanation  might  be 
radiative  damping.  A conclusive  explanation  for  the  large 
departure  of  over  ocean  temperature  spectra  from  C - ^ is 

still  forthcoming,  but  seems  to  be  associated  with  a co 
spike"  signature  in  the  temperature  signal  at  the  interface  of 
a large  scale  roll-like  structure  over  the  open  ocean  during 
unstable  conditions.  "Cold  spikes"  appear  to  be  due  to  some 
sort  of  evaporative  phenomenon  - but  the  precise  mechanism  is 
not  understood  at  this  time,  and  is  under  active  investigation. 
Simultaneous  measurements  of  the  water  vapor  concentration 
showed  no  anomalous  behavior  as  found  for  the  temperature.  ^ 
Further  discussion  of  this  phenomena  is  included  in  Dreyer  s 
dissertation  and  in  Friehe  et  al  (attached)  . 
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Turbulence  Measurements  in  tie  Marine  Boundary  Layer 
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1 . Introduction 

As  approximately  three-quarters  of  the  earth's  surface 
is  covered  by  water,  about  three-quarters  of  the  atmospheric 
boundary  is  over  a water  surface.  The  transfer  rates  of 
momentum,  heat  and  moisture  through  the  marine  bcundarv  laye. 
have  large  effects  on  the  global  weather  system,  oceanic  and 
atmospheric  circulations  and  heat  contents  and  their  inter- 
actions. However,  more  effort  has  been  expended  in  the  past 
on  atmospheric  measurements  over  land  than  over  the  seas. 
Recently  (since  about  1968)  extensive  and  detailed  measurements 
have  been  made  in  the  marine  boundary  layer.  it  is  our 
purpose  here  to  describe  the  methods  used  and  measurements 
obtained  by  our  research  efforts  in  the  marine  boundary  layer 
in  the  last  few  years. 

The  marine  boundary  layer  is  char ac ter ized  by  turbulent, 
rather  than  laminar,  motion.  Therefore,  in  making  measure- 
ments in  the  marine  boundary  layer,  one  is  making  measurements 
of  the  turbulent  velocity  field  and  associated  fields  of 
temperature  and  water  vapor  content  that  are  mixed  by  the 
turbulent  motion.  In  order  to  measure  velocity,  temperature 
and  water  vapor,  instrumentation  of  high  sensitivity  and 
bandwidth  is  required  to  measure  the  entire  spectrum  of  the 
turbulent  signals.  For  example,  the  energy  of  the  vertical 
component  of  the  turbulent  velocity  vector  is  typically 
concentrated  at  about  0.01  to  0.1  Hz,  and  is  dissipated  by 
viscosity  at  between  100  to  1000  Hz  for  the  range  of  usual 
marine  boundary  layer  conditions.  Thus  to  measure  the  entire 
bandwidth  of  the  signal  requires  sensor  (and  recording  and 
associated  instrumentation)  bandwidth  from  essentially  d.c. 
to  about  1 kHz.  Signal  to  noise  requirements  are  also  large, 
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since  the  power  spectra  of  most  of  the  signals  have  a fairly 
/apid  roll-off;  in  the  inertial  subrange  power  spectrum  is 
proportional  to  frequency  to  the  -5/3  power.  The  above 
remarks  apply  to  measurements  over  land  as  well  as  over  the 
oceans.  Ocean  measurements  are  further  complicated  by 
additional  factors,  such  as  typically  lower  signal  levels,  a 
more  s?vere  environment,  and  difficulty  in  obtaining  measure- 
ments that  are  not  affected  by  wave-induced  motion  of  the  ship 
or  instrument  platform.  Also,  over  land  the  measurement  of 
moisture  is  not  often  made,  or  measurements  are  made  in  areas 
or  under  conditions  under  which  moisture  is  not  important. 

Over  the  oceans,  moisture  is  obviously  important,  but  a 
difficult  variable  to  measure. 

In  summary,  measurements  in  the  marine  boundary  layer 
provide  information  about  the  exchange  of  heat,  moisture  and 
momentum  between  the  oceans  and  atmosphere.  These,  together 
with  radiation  and  precipitation  measurements,  will  further 
our  understanding  of  weather,  oceanic  processes  and  their 
interactions.  Also,  the  marine  boundary  layer  is  a "laboratory' 
for  fundamental  turbulence  measurements.  Most  turbulence 
theories  apply  at  large  Reynolds  numuers  (the  Reynolds  number 
is  the  ratio  of  inertial  to  viscous  forces  in  fluid  motion), 
and  in  the  atmospheric  boundary  layer  the  Reynolds  numbers 
are  typically  several  orders  of  magnitude  greater  than  that 
obtainable  in  laboratory  flows.  The  scales  of  the  notion 
are  also  large  so  that  sensor  sizes  are  at  least  feasible  to 
manufacture  to  resolve  the  smallest  scales  in  the  flow  field. 


2.  Measurement  Facility:  FLIP 

Measurements  of  atmospheric  turbulence  over  the  open 
ocean  have  been  made  from  at  least  four  types  ot  platforms: 
moving  ships  (Mitsuta,  et  al  (1974)),  a gy r o - s t ab il iz ed  buoy 
(Dunckel , et  al  (1974)),  a aircraft  (Miyake,  et  al  (1970)), 
and  from  the  manned  spar  buoy  R/P  FLIP*  (Pond,  et  al  (1971), 
Gibson  et  al  (1970),  for  example).  In  this  article  we  will 
review  primarily  measurements  made  from  FLIP.  FLIP,  as 
configured  for  a recent  atmospheric  experiment,  is  shown  in 
Figure  1.  About  100  m of  FLIP  is  underwater  so  that  surface 
wave  generated  pressures  have  little  effect  on  FLIP'S  vertical 
motion.  (For  more  detailed  information  on  the  design  and 
performance  of  FLIP,  see  Rudnick  (1967)).  Atmospheric 
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sensors  are  placed  at  the  end  of  the  large  port  boom  to  avoid, 
as  much  as  possible,  effects  due  to  the  distortion  of  tie  wind 
field  around  FLIP.  The  starboard  boom  has  a sail  on  it  to 
balance  the  drag  of  the  larger  boom  so  that  FLIP's  keel  is 
kept  pointing  into  the  wind.  (In  Figure  1,  the  wind  would  be 
coming  out  of  the  picture  toward  the  reader.)  Sensors  are 
located  on  small  masts  above  boom  level  , or  below  to  obtain 
measurements  near  the  surface  (not  shown  in  Figure  1).  Some 
measurements  at  about  30  m above  the  mean  ocean  surface  have 
also  been  made  from  a collapsible  vertical  mast  (shown  lowered 
in  Figure  1).  Signal  cables  are  strung  from  the  sensors  to 
the  laboratory,  located  at  i he  level  of  the  booms.  About  12 
"bays"  of  electronic  instrumentation  can  be  placed  in  the 
laboratory,  and  about  10  scientists  can  be  accommodated  for  a 
maximum  single  period  of  30  days. 

The  motion  oi  FLIP  is  of  course  much  less  than  the  motion 
of  a regular  surface  ship.  However,  the  residual  motion  does 
contaminate  some  measurements,  especially  that  of  the  vertical 
velocity  component  . Our  method  of  correcting  for  this  mot  on 
is  described  later. 


3.  Experimental  Methods 
3.1  Instrumentation 


Velocity 


Three  independent  devices  are  generally  used 
to  measure  the  wind  velocity.  A standard  micrumeteor ological 
cup  anemometer  is  used  to  record  the  horizontal  wind  speed  and 
fluctuations  up  to  about  0.5  Hz.  A three-dimensional  somc^ 
anemometer  is  used  to  meat  ire  simultaneously  the  three  ccmpo 
Tents  of  the  velocity  vector  from  d.c.  to  10  Hr.  The  cup  and 
sonic  anemometers  do  not  require  calibration  in  the  field,  but 
are  tested  in  the  UCSD  AMES  low  speed  wind  tunnel  to  verity 
operation  before  and  after  field  trips.  For  high  frequency 
velocity  component  measurements,  linearized  constant  tempera 
anemometers  are  used  with  hot-wire  or  hot-filmsingle  or  double 
sensor  probes.  The  double  sensor  probes  are  in  an  X con 
figuration  for  the  measurement  of  the  horizontal  and  vertical 
velocity  components.  Bandwidth  is  typically  d.c.  to  several 
kHz  with  probe  cables  50  m long.  The  probes  are  calibrated 
in  a small  variable  speed,  controlled  air  temperature  wind 
tunnel  inside  the  laboratory  on  FLIP  before  and  atter  use, 
without  disconnecting  the  probes  from  the  anemometer  circuitrv. 


<°-°25  "1dU",'r.;r  ih.n  h«-5iti  sensor s (0.0038  mm 
r„e  ..rlne  .tm„,pK«r«  b.  ter  th.nahot  (1„*)>.  The 

diameter,  1.2  mm  lengtl J equipped  with  a protractor 

Cao“bd:^ne  -»»»“<>"  o,  Che 

X probes . 


b.  Temperature 
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c . Water  Vapor 


Measurement  of  the  fluctuations  of  water  vapor 
content  with  spatial  and  temporal  resolution  equal  to  that 
for  velocity  and  temperature  fluctuations  is  not  possible 
with  existing  instrumentation.  Wetted  thermocouples  or 
thermistors  have  slow  response  times,  and  the  p s y c h r ome t r ic 
relationship  has  to  be  used  to  obtain  measurements  of  water 
vapor  density  (absolute  humidity,  pgm/cm3).  We  have  chosen 
to  use  the  Lyman-alpha  humid  i ome  to  r , which  directly  measurer, 
ter  vapor  density  in  the  atmosphere.  (See  A.  Buck  (1973), 
for  a detailed  description  of  this  instrument.)  1 he 

instrument  measures  the  absorption  of  the  Lyman-alpha  line  of 
hydrogen,  which  is  assumed  to  be  caused  only  by  the  hydrogen 
contained  in  wate’1  vapor  in  the  measuring  path.  The  absorption 
approximately  follows  Beer's  law,  so  that  on  exponential 
dependence  of  detector  tube  current  on  water  vapor  concentration 
is  obtained.  The  output  is  often  linearized  with  a logarith- 
mic amplifier.  The  commercial  instrument  has  hydrogen  source 
tubes  and  nitrous  oxide  detector  tubes  of  about  2 cm  in  diameter 
and  5-8  cm  long.  The  path  length  is  adjustable  from  about  0.5 

cm  to  5 cm.  Thus  the  spatial  resolution  is  not  adequate  to 
measure  the  "dissipative"  scales  of  the  fluctuating  water  vap^r 
field,  of  the  order  of  1 mm.  A further  operational  problem 
with  the  Lyman-alpha  device  is  that  the  tube  windows  are  of 
magnesium  flouride,  which  transmit  the  Lyman-alpha  line 
1216  A),  but  are  slightly  soluble  in  water.  Attenuation  of 
the  signal  occurs  due  to  deterioration  of  the  windows  from 
spray  or  exposure  during  squalls.  Unfortunately,  a better 
window  material  has  not  been  found. 

To  measure  average  and  low  frequency  water  vapor 
content  and  to  calibrate  in  situ  the  Lyman-alpha  humid iometer , 
a "cooled  mirror"  dew  point  instrument  is  used.  This  device 
has  a dew  point  accuracy  of  about  0.3  C (about  0.2  vigm/cm3) 

The  basic  principle  of  operation  is  the  measurement  of  the 
temperature  of  a metal  membrane  that  is  heated  or  cooled  until 
dew  is  formed  on  the  surface. 


d . Motion 

Slight  motion  of  the  FLIP  affects  primarily  the 
measurement  or  th  instantaneous  vertical  velocity  component. 
The  outputs  of  the  sonic  and  X-probe  anemometers  are  the  same 
for  an  actual  variation  of  the  turbulent  velocity  vector  or  an 
equivalent  motion  of  the  sensors  in  a uniform  flow.  Since 
the  horizontal  velocity  component  is  typically  an  order  of 
magnitude  greater  than  the  vertical  component,  small  tilting 
angles  about  the  vertical  axes  of  the  sensors  can  cause 
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appreciable  contamination  of  the  vertical  velocity  component 
by  the  horizontal  component.  Also,  there  are  true  velocit> 
effects  in  the  measured  signals  due  to  the  fact  that  iLIP  is 
usually  slowly  moving  with  respect  to  a fixed  reference 
frame.  To  correct  the  first  effect  of  tilt  sensitivity  of 
the  anemometer  sensors,  a vertical  gyroscope  (originally 
obtained  by  R.  Davis  of  SIO)  is  placed  next  to  the  sensors 
to  measure  the  instantaneous  tilt  angle.  The  tilt  correction 
is  performed  on  the  data  by  digital  analysis  in  the  (land) 
laboratory.  The  total  motion  of  FLIP  is  in  general  compli- 
cated, since  there  a.c  three  components  of  angular  motion 
and  three  of  translational  motion.  An  airline  mertiai 
navigation  unit  (very  similar  to  the  one  used  on  the  NCAR 
Electra)  has  been  used  once  to  measure  the  six  components, 
and  the  data  is  presently  being  reduced. 


3.2  Instrumentation  Package 

The  total  instrument  package  is  shown  in  Figure  2, 
with  the  exception  of  the  dew  point/mean  air  temperature 
unit.  The  mounting  frame  is  quite  rigid  and  the  angles  be- 
tween  the  gyroscope  and  anemometer  sensor  mounting  pads  are 
precisely  set  in  the  laboratory.  The  gyroscope  is  calibrated 
on  a special  tilt  table  in  the  (land)  laboratory  with  respect 
to  the  vertical. 


3.3  Signal  Processing  and  Recording 

Data  are  recorded  on  FLIP  mainly  on  analog  FM 
instrumentation  magnetic  tape  recorders.  Many  channels  (up 
to  20),  several  of  high  bandwidth,  must  be  recorded  simulta- 
neously. The  information  rate  is  too  large  for  practica 
digital  recording,  although  in  some  experiments  we  have  used 
a on-line  minicomputer  system  wliich  stores  data  on  digital 

tape  . 

Since  the  root-mean-square  signal  to  noise  ratio 
of  FM  recorders  is  only  about  50  dB,  pre-processing  of  the 
basic  signals  is  required  before  recording  to  preserve  their 
full  dynamic  ranges,  which  may  exceed  50  dB.  S ’ gnal 
processing  equipment  consists  of  precision  d.c.  offset  and 
gain  amplifiers,  analog  filters,  differentiators  and  pre- 
emphasis  circuits  (of  a specially  shaped  transfer  function  to 
r re-whit  e n typical  atmospheric  turbulent  velocity  ard 
temperature  spectra).  Time  division  multiplexing  is  used  lO 
record  low  frequency  information  from  several  (up  to  /) 
instruments  on  one  analog  tape  track.  Observationa 
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equipment  used  on  board  FLIP  to  monitor  signal  quality 
includes  usual  storage  oscilloscopes  anti  strip  chart 
recorders,  and  a power  spectrum  analyzer.  As  mentioned 
above,  a minicomputer  has  been  used  to  collect  data,  and 
also  perform  in  situ  data  checks. 


3.4  Digital  Analysis 

The  recorded  data  is  converted  to  digital  form  in 
the  laboratory  at  UCSD  (Gibson  (i973)).  Two  computer 
systems  are  used  for  various  statistical  and  spectral 
analyses;  IBM  1130  and  TI  980  in  the  laboratory  or  the 
University  CDC  3600.  Digitization  for  the  1130  and  980  are 
14  bit  with  sampling  rates  up  to  8 kHz  per  channel;  for  the 
3600,  12  bits  are  used.  Four  channels  of  matched  48  dB/ 
octave  analog  low  pass  filters  are  available  to  minimize 
aliasing  for  spectral  analysis.  Spectral  analysis  is  done 
using  fast  Fourier  transform  routines. 


4.  Results  and  Discussion 

A typical  time  series  trace  of  horizontal  and  vertical 
components  of  the  velocity  vector,  temperature  and  humidity 
fluctuations  is  shown  in  Figure  3.  The  data  were  obtainei 
from  FLIP  at  a height  of  about  10  m over  the  water  under 
condit.ions  of  relatively  cold,  dry  air  over  hot  (wet)  water. 
The  fluxes  of  momentum,  heat  and  water  vapor  can  be  directly 
calculated  from  the  time  series.  The  vertical  flux  of 
momentum  is  proportional  to  the  covariance  of  the  horizontal 
and  vertical  velocity  components;  the  vertical  flux  of 
sensible  heat  is  proportional  to  the  covariance  of  the 
fluctuating  temperature  and  vertical  velocity;  and  the 
vertical  flux  of  moisture  is  proportional  to  the  covariance 
of  the  fluctuating  water  vapor  concentration  and  the  vertical 
velocity.  The  evaporation  of  water  at  the  air-sea  interface 
represents  a latent  heat  loss  from  the  ocean.  For  the 
conditions  described  above,  the  covariance  of  the  two  velocity 
components  is  negative,  while  the  covariances  of  the  vertical 
velocity  component  with  the  scalar  fields  are  positive.  The 
signs  of  the  covariances  are  somewhat  difficult  to  determine 
by  inspection  of  the  time  series  plots  since  the  correlation 
coefficients  (absolute  value)  are  usually  less  than  one  — half. 
There  is  a fairly  high  correlation  between  the  temperature 
and  humidity  fluctuations,  and  the  time  trace  exhibits  a 
characteristic  ramp  structure  of  slowly  increasing  temperature 
and  humidity  with  time,  leading  to  a sharp  interface.  The 
correlation  of  temperature  and  humidity  is  not  always  perfect: 
at  some  of  the  ramp  interfaces  the  temperature  decreases  to 


abnormally  low  values  ("cold  spikes")  compared  to  the  behavior 
of  the  humidity.  An  example  of  such  a time  series  is  shown 
in  Figure  4.  The  anomalous  "cold  spikes"  of  the  temperature 
field  were  apparently  first  observed  in  project  BOMEX  (Pond, 
et  al  (1971)).  Hypotheses  advanced  to  explain  the  "cold 
spikes"  include  radiative  effects  and  mixing  of  upper-level 
air  of  different  temperature  and  humidity  profiles  (Pond  et 
al  (1971)),  and  evaporation  of  spray  at  the  ramp  interfaces 
(Dreyer  (1974)).  There  have  not  been  reports  of  "cold 
spikes"  in  temperature  fields  observed  over  land. 

Spectral  analysis  has  been  of  great  importance  in 
understanding  the  physics  of  turbulent  motion  and  mixing  of 
scalar  fields  by  turbulence.  The  similarity  hypotheses  and 
inertial  subrange  prediction  of  Kolmogorov  should  be  partic- 
ularly applicable  to  the  high  Reynolds  number  turbulent 
atmospheric  boundary  layer.  Power  spectra  of  velocity 
components,  temperature  fluctuations  and  humidity  fluctuations 
from  low  (*10-3Hz)  to  moderate  ('10  Hz)  frequencies  obtained 
from  FLIP  are  shown  in  Figures  5a  and  5b.  The  velocity  data 
were  obtained  with  a sonic  anemometer,  and  show  an  approach 
with  increasing  frequency  to  the  Kolmogorov  inertial  subrange 
form  of  spectral  energy  proportional  to  the  -5/3  power  of 
frequency.  (The  Kolmogorov  inertial  subrange  form  is 
expressed  in  wave  number;  the  conversion  from  frequency  is 
made  by  equating  the  streamwise  spatial  coordinate  to  the 
product  of  time  and  the  mean  velocity  past  the  fixed  measuring 
probe.)  The  humidity  fluctuation  data  obtained  with  the 
Lyman-alpha  device  also  shows  an  inertial  subrange  behavior 
for  scalars,  as  predicted  by  the  Obukho v-Cor r s in  analyses 
similar  to  Kolmogorov's  for  the  velocity  field.  However, 
the  temperature  fluctuations  do  not  exhibit  an  inertial  sub- 
range. The  slopes  of  the  temperature  power  spectra  for  the 
marine  boundary  layer  are  consistently  less  than  -5/3  over  a 
wide  frequency  range,  apparently  largely  due  to  the  cold 
spikes".  The  basic  reason  (or  reasons)  for  this  behavior 
are  not  fully  understood.  To  a first  approximation,  tem- 
perature and  water  vapor  content  should  behave  as  identical 
passive  scalars.  The  molecular  diffusivity  of  heat  in  air. 
is  within  20%  of  that  of  water  vapor  in  air,  and  the  governing 
equations  are  the  same,  except  for  source  or  sink  terms  due 
to  radiative  effects  and  evaporation  of  spray. 

High  frequency,  fine  scale  results  for  the  velocity  and 
temperature  fields  are  presented  as  power  spectra  of  the  time 
derivatives  of  the  variables.  The  frequency  dependence  in 
the  inertial  subrange  for  the  power  spectrum  of  a derivative 
becomes  f-5/3*  f2=fl/3,  where  f=frequency.  As  shown  in 

Figure  6,  the  velocity  power  spectrum  exhibits  the  f form, 

whereas  the  temperature  derivative  obtained  at  three  heights 


does  not,  except  for  a small  region  of  tangency  before  the 
diffusive  roll-off.  To  confirm  the  departure  of  the  over 
ocean  temperature  results  from  the  inertial  subrange  form, 
derivative  power  spectra  were  obtained  in  a laboratory  flow 
at  moderately  high  Reynolds  number  - a heated  turbulent  jet 
(so  large  that  the  experiment  had  to  be  performed  in  the 
UCSD  gymnasium).  Essentially  the  same  equipment  for  the 
field  experiments  was  used,  and  both  the  temperature  and 
velocity  derivative  spectra  substantially  followed  the 
inertial  subrange  form  indicating  that  anrmalous  instrument 
response  was  not  a factor.  It  has  also  been  determined  that 
the  fine  scale  structure  is  not  isotropic,  as  evidenced  by 
non-zero  values  of  the  skewness  of  the  time  derivatives  of 
temperature  in  the  marine  boundary  layer  as  found  by  Gibson, 
et  al  (1970)  and  Boston  and  Burling  (1972).  The  universal 
constant  associated  with  the  dimensionless  form  of  the 
inertial  subrange  formula  for  temperature  was  found  in  those 
two  studies  to  be  greater  (by  factors  of  2 to  3)  than  values 
measured  in  previous  laboratory  studies. 


5 . 


Conclusions 


Turbulence  measurements  in  the  marine  boundary  layer 
provide  insight  into  the  physical  mechanisms  of  air-sea 
interaction  and  high  Reynolds  number  turbulence.  Measurements 
of  the  fluxes  of  momentum,  sensible  heat  and  water  vapor  over 
the  oceans  are  required  to  further  our  understanding  and 
modeling  of  atmospheric  and  oceanic  processes.  Direct  flux 
measurements  (calculating  the  appropriate  covariances  of  the 
turbulent  signals)  are  essential  to  confirm  or  calibrate  other 
indirect,  simpler  techniques.  Such  efforts  are  convenient  y 
performed  from  a stable  platform  such  as  the  1LIP.  The 
simplest  indirect  flux  estimate  method  is  the  use  of  the  bulk 
aerodynamic  formulas  (Roll  (1964))  which  relate  the  fluxes  to 
observations  of  mean  air  and  sea  temperatures,  mean  water  vapor 
concentration,  and  mean  wind  speed.  The  empirical  coef  l- 
cients  in  the  formulas  are  determined  experimentally  by 
compaLison  to  simultaneous  direct  flux  measurements  as  by  Pond 
et  al  (1971),  for  example.  Another  technique  involves  fine 
scale  measurements  of  the  dissipation  rates  of  velocity  and 
temperature  fluctuations  (the  dissipation  rates  are  proportional 
to  the  variances  of  the  time  derivatives  of  the  fluctuations). 
From  the  measured  rates,  the  momentum  and  sensible  heat  fluxes 
may  be  inferred  (Stegen,  et  al  (1973)).  The  simpler  tech- 
niques are  of  interest  as  they  do  not  involve  measurement  of 
the  vertical  velocity  component,  and  therefore  may  work  on 
conventional  ships. 
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Fine  scale  measurements  in  the  marine  boundary  layer 
have  revealed  some  departures  from  theoretical  predictions 
and  assumptions  about  turbulence  at  high  Reynolds  number, 
particularly  for  the  temperature  field.  The  reasons  for 
the  departures  are  not  presently  understood  and  are  an 
active  subject  of  present  research. 
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Figure  la  and  lb: 


FLIP,  "flipped"  and  being  rigged 


for  marine  boundary  layer  measurements.  FLIP  is  towed, 
in  the  horizontal  position,  to  and  from  the  experimental 
site  as  in  Figure  lb.  FLIP  does  not  have  any  propulsion 
power  of  her  own,  is  properly  called  a spar  buoy.  Transition 
from  the  horizontal  to  the  vertical  takes  about  30  minutes; 
from  vertical  to  horizontal  about  15  minutes.  The  side  booms 
are  folded  for  the  horizontal  position.  Electronic 
instrumentation  is  mounted  sideways  in  standard  racks  in 
the  laboratory  in  the  horizontal  position.  Transfer  of 
scientists  and  small  equipment  to  support  vessels  is  possible 

by  use  of  the  small  boat. 
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Figure  2a:  Turbulence  and  micrometeorological 


instrument  package.  A.  Velocity  sensors:  (1)  Three- 


dimensional sonic  anemometer,  (2)  Hot-wire/ f ilm  velocity 


probe  (single  or  X-sensors)  for  use  with  constant  temperature 


anemometers,  (3)  Cup  anemometer.  B.  Temperature  sensors: 


(4)  Platinum  wire  temperature  sensor  probe  for  use  with  the 


a.c.  bridge  circuit,  (5)  Thermistor  probe.  C.  Water  vapor 


sensors:  (6)  Lyman-Alpha  humid iometer . (Not  shown  is  mean 


air  temperature  and  dew  point  instrument).  D.  Motion  sensor: 


(7)  Vertical  gyroscope. 


Figure  2b:  Front  view. 
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Figure  3:  Time  series  of  horizontal  and  vertical 

velocity  components,  temperature  and  humidity.  The  vertical 
flux  of  sensible  heat  is  proportional  to  the  covariance  of  the 
temperature  and  vertical  velocity  time  series;  the  vertical 
flux  of  moisture  is  the  covariance  of  the  absolute  humidity 
and  the  vertical  velocity;  the  vertical  flux  of  momentum  is 
proportional  to  the  covariance  of  the  horizontal  and  vertical 
velocity  components. 
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Figure  5a:  Power  spectra  of  turbulent  velocity 

fluctuations  obtained  over  the  open  ocean  with  a three- 
dimensional  sonic  anemometer.  Slope  of  the  solid  line  is 
-5/3  for  the  inertial  subrange.  The  slight  peaks  in  the 
spectra  at  “10”^  Hz  may  be  caused  by  surface  waves. 


Figure  5b:  Power  spectra  of  temperature  and  humidity 

fluctuations  obtained  with  a platinum  cold  wire  sensor  and 
a Lyman-a lpha  humid iome t er . The  humidity  power  spectrum 
follows  the  -5/3  slope,  but  the  temperature  spectrum  does 
not.  Overlapping  spectral  bands  (two  for  humidity  and  three 
for  temperature)  were  used  to  cover  the  large  dynamic  ranges 
of  the  signals.  The  frequency  response  of  the  temperature 


sensor 


is  greater  than  that  of  the  humidity  sensor. 
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Figure  6:  Spectra  of  the  time  derivatives  of  velocity 

and  temperature.  The  data  at  3m  height  were  obtained  with 
a hot  wire  (velocity)  -Pt  cold  wire  (temperature)  probe  with 
1 mm  sensor  separation.  The  Kolmogorov  inertial  subrange 
is  given  by  a slope  of  +1/3  for  the  power  spectrum  of  a 
derivative.  The  velocity  exhibits  an  inertial  subrange, 
whereas  the  temperature  does  not. 
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„ For  a passive  scalar  at  high  Reynolds  number  and  at  a 

Prandtl  number  (Pr  = -)  near  unity  the  temperature  derivative 

D 

spectrum  should  have  a universal  inertial  subrange  of  slope 

+ 1/3.  The  predicted  form  for  the  derivative  spectrum  is 

$ • (k)  = &xe  i/3^1/3  where  x is  the  scalar  dissipation  rate, 

T 

e the  viscous  dissipation  rate,  k is  the  wavenumber  and  6 a 
"universal"  constant.  In  the  atmospheric  boundary  layer  B 
is  not  constant  and  over  the  open  ocean  the  slope  of  the 
derivative  spectrum  varies  between  0.7  and  the  predicted  value 
of  0.33.  These  variations  in  slope  seem  to  be  related  to  the 
background  meteorological  conditions,  steeper  slopes  occurring 
under  cool,  dry  conditions  and  slopes  nearer  =1/3  under  warm, 
humid  conditions.  The  evaporation  of  spray  in  the  near 
water  zone  could  cause  variations  in  slope  in  conjunction 
with  radiational  effects  due  to  high  absolute  humidities  over 
the  ocean. 
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